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Abstract 

The desirability of evaluating the effective potential in field theories near a 
phase transition has been recognized in a number of different areas. We show 
that recent Monte Carlo simulations for the probability distribution for the 
order parameter in an equilibrium Ising system, when combined with low- 
order renormalization group results for an ordinary (f/^ system, can be used to 

extract the effective potential. All scaling features are included in the process. 
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I. INTRODUCTION 



In the realm of statistical physics, as well as in quantum field theory, the need for non- 
perturbative approaches to study the vicinity of a phase transition has long been appreciated. 
Results of recent simulations [1-4] show that the form of the effective potential for a variety 
of systems can be determined using Monte Carlo methods. Such methods have also been 
applied to the study of the electroweak phase transition via lattice Monte Carlo methods [5] . 
Additional continuum limit work has recently been performed on a X(f)^ theory in three- 
dimensions [6]. Extracting the effective potential provides one with an additional window 
onto the equilibrium - and possibly the dynamical - behavior of a system near and at 
its critical point. One use of Monte Carlo methods is to test scaling forms based on the 
Renormalization Group. A second is to explore the equilibrium behavior associated with 
such phenomena as two-phase coexistence and broken continuous symmetry. 

Binder [7] and more recently Tsypin [1] and Chen and Dohm [8] have utilized Monte 
Carlo and binning techniques to generate probability distributions of the order parameter 
of an 0(1) (Ising-like) system in the vicinity of its critical point. In Tsypin's work data 
were generated at a single temperature, and over a range of symmetry-breaking fields in 
order to explore the full scaling domain. Above the critical temperature the distribution 
was found to be consistent with a low-order polynomial form for the effective potential, 
which, as discussed below, is essentially the logarithm of the probability distribution. The 
results of simulations were most closely fit when the effective potential was terminated at 
the sixth order. The standard quartic Ginzburg-Landau form for the effective potential 
failed to produce an adequate fit for any choice of coefficients, and the addition of eighth 
and higher order terms did not materially improve the agreement between the fitting form 
and the results of simulations. 

This remarkable set of results persists below the critical temperature [2]. Here, it was 
found that a sixth-order polynomial form for the effective potential leads to an outstanding 
match to the results of simulations. In this case the best-fit effective potential was found to 



2 



contain a negative quartic term. An additional feature of Tsypin's fitting form is a prefactor, 
proportional to the square root of the second derivative of the effective potential. 

It is reasonable to ask whether or not the data are consistent with a scaling form for the 
effective potential. In fact, the nature of the polynomial strongly suggests a form incorpo- 
rating thermodynamic scaling with exponents appropriate to the Ising model. The rationale 
for this will be explored in more depth in subsequent sections. 

It should be noted that the renormalization group has been utilized previously to generate 
order parameter distributions in the vicinity of the critical point [4,8,9]. The principal 
novelty in this paper is the attempt to produce a unified form for the order parameter 
distribution that fits data above, below and at the critical point. It is found that the 
inclusion of the prefactor mentioned above leads for excellent agreement with data based on 
Monte Carlo simulations. The prefactor does not play an important role in the fitting of a 
scaling form to the order parameter distribution in the high and low temperature phases. 
However, this finite size correction is essential to the construction of a high-precision fit to 
data at the critical point. 

Checks on the quality of the fit include comparisons of a universal ratio of moments of the 
distribution with values in the literature determined by a variety of alternative techniques, 
and comparison between our optimum fit value of the fourth order coupling u with a variety 
of other determinations of that quantity. These checks are very encouraging. It appears 
that a scaling form based on renormalized mean field theory agrees to a high degree of 
accuracy with published results. On this basis, one has confidence that an equation of state 
constructed in a similar manner from renormalized mean field theory will accurately describe 
systems in the or Ising model, universality class. 

The remainder of this paper is laid out as follows. In the next section the effective 
potential is defined and a phenomenological scaling description is presented. In Sec. Ill a 
renormalization group derivation of the scaling form is given, and in the following section 
the "prefactor" is discussed. In Sec. V finite size effects are addressed, while Sec. VI assesses 
the success of the scahng form of the effective potential as a fit to Tsypin's Monte Carlo 
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simulations. Sec. VII is devoted to concluding remarks. An Appendix contains a discussion 
of the effect of the prefactor on an important universal quantity. 



II. SCALING CONSIDERATIONS 

We imagine a scalar (Ising-like) system described by order parameter, and ordering 
field, h, in equilibrium at temperature T. If the system is described by Hamiltonian 
the constrained (reduced) free energy, F{T, h; M), is given by 

Z = exp[-F{t,h;M)]= J D(j) 6 {M - I{(j)))exp[-H{(j))/T] , (2.1) 

where /(0) = J (j) df^x is the average order in configuration {0}, represents the volume 
of the system, and d is the dimensionality. The effective potential Ti, is usually defined as 

Z = exp(-H) (2.2) 

which identifies it with the constrained free energy. According to standard scaling notions, 
when the linear extent, L, of the system is sufficiently great, the effective potential, 7Y, of 
an 0(1) system will have the form 

n {t, h, M) = \M\^'i/i<i-''+^) T (t |M|-2/-(<^-2+'?) , h |^|-(d+2-r,)/{<i-2+^)>j _ ^2.3) 

The parameter t is the reduced temperature, t — {T — Tc)/Tc, where Tc is the critical tem- 
perature, and the ordering, or symmetry-breaking, field h has been introduced above. The 
quantity L is a characteristic linear extent of the system, and d is its spatial dimensional- 
ity. The exponents v and r] appear in the (unconstrained) two-point correlation function 
C(f-P,t,/i), 

C {r, t, h) = 1^1"^'"'^'^ C (\r\ t\ hr<'+''-^^'^) . (2.4) 

All thermodynamic exponents follow from the two correlation function exponents, u and r^, 
given the standard scaling and hyperscaling relations [10]. When the spatial dimensionality 
is d = 3, the critical exponents are [11] 



1/ = 0.63 ± 0.002 (2.5) 
= 0.037 ±0.001 (2.6) 

An immediate consequence of expression (2.3) is that the effective potential at the critical 
point {t = h = 0) has the form 

^^critical oc L'^ \Mf'^/^'^~^+''^ = L'^ |M|^'-^ . (2.7) 

It is reasonable to expect that, within the scaling regime, at low enough reduced temperature 
and ordering field, and for appropriate values of the order parameter (that is, if the order 
parameter controls the decay of correlations) the effective potential will behave in the same 
way. If one attempts to approximate the effective potential by a polynomial in the order 
parameter, M, the "best" one will, most likely, not have terms that go beyond sixth order. 
Interestingly, this is precisely the result of Tsypin's unbiased attempt to fit his simulation 
data to a polynomial effective potential: A fourth order, Ginzburg- Landau form fails to 
adequately represent the data, while the coefficients of terms beyond sixth order are so small 
as to cast doubt on their appearance in the true effective potential. This result holds both 
above and below the critical temperature — that is, whether the system is in the disordered 
or the ordered state. 

Now, the form (2.3) is expected to apply only when the order parameter M lies at - 
or close to - the value which minimizes the free energy at fixed t,h. In particular, there 
are known difficulties in applying the above scaling form (2.3) to the case of a system 
in the coexistence region, which impinge directly on the scaling regime. The equilibrium 
state within the coexistence regime will generally involve one or more interfacial regions 
separating two homogeneous thermodynamic phases. Because of this, the scaling form, 
which is hypothesized for a single homogeneous phase, does not necessarily do an adequate 
job of describing the behavior of a system in which two different phases coexist. On the other 
hand, a scaling form ought to predict with great accuracy the dependence of thermodynamic 
functions on temperature and ordering field as these fields approach the coexistence line, 
that is, at the phase boundary for two-phase coexistence. 
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There is one more consequence of scaling that is worthy of note. The scahng form in 
Eq. (2.3) together with correlation function scaling implies the following form for the free 
energy [12,13] 

= G (|M|t|-^('^-2+'')/2, Lltl", ^ j . (2.8) 

In the first line of Eq. (2.8) the quantity ^ is the bulk correlation length a,t h — 0. 

III. RENORMALIZATION-GROUP BASED FORM FOR THE EFFECTIVE 

POTENTIAL 

The general form in Eq. (2.3) restricts, but does not specify, the detailed dependence 
of the effective potential Ti. on t and h. Other considerations are required for an explicit 
evaluation. One candidate form is based on the Renormalization Group, particularly the 
field-theoretical expansion in e = 4 — d, where, as above, d is the spatial dimensionality of 
the system. We consider a standard cf)^ Hamiltonian 

^-Jll + r<l>') + u<l>'] d'x . (3.1) 

At lowest non-trivial order in e, the effective potential has the form [14] 



n{M,t,h)^e 



-l*d 
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(3.2) 



This approximate effective potential is at the level of a renormalized Ginzburg-Landau free 
energy, and is hence equivalent to "renormalized mean-field theory" . The renormalization, 
indeed the full effect of the Renormalization Group trajectories, resides in the quantity e^* 
which is equivalent to the "block spin" size. This key quantity may be determined via 

The effective potential described by Eqs. (3.2) and (3.3) is fully consistent with the scal- 
ing hypotheses embodied in Eq. (2.3). It is is completely determined once the coupling 
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constant u has been set, along with "metrical" factors associated with the scales of h and 
t. The inclusion of higher order, loop expansion contributions leads to modifications of the 
effective potential, so, on the face of it, there is no reason to expect that (3.2) will be exact. 
However, given that it is entirely consistent with critical point scaling, that it embodies the 
full RG trajectories, and that it represents the leading order contribution in a systematic 
Renormalization Group expansion for the free energy, it seems at least a plausible "zeroth 
order" candidate for the critical effective potential. This optimism is bolstered by other 
situations in which "renormalized mean-field theory" provides surprisingly good numerical 
results when applied in three dimensions. [15] 



IV. THE PREFACTOR 

In this section we discuss the prefactor introduced by Binder [16] and used by Tsypin in 
fitting Ising simulations [1,2]. A variety of arguments have been advanced for the existence 
of this term. Here we introduce additional justification for the prefactor. 

Suppose we allow spatial variations of the order parameter around the globally con- 
strained value M in the form M + a{k). Note that {cr{k)) = 0. Neglecting terms of greater 
than quadratic order in the a{kys, the dependence of the effective Hamiltonian on the a{kys 
at gaussian order will be 



a{k) 



T 



1 



1 



E 2 + ^n"a{k)a{-k) 
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(4.1) 



The term H" is the second derivative with respect to the order parameter M of the effective 
potential Ti. at renormalized mean field level discussed above. As noted the sum excludes 
the A; = mode of the order parameter fluctuations. Integration over the cr(A;)'s yields the 
following "gaussian" correction to the effective potential 

As above VL is the spatial volume of the system. The last term on the right hand side of Eq. 
(4.2) is the leading order difference between the sum on the left hand side and its asymptotic 



limit as the integral on the right. This term can be thought of as a finite size correction to 
the extensive result in the limit, Q ^ oo. A general feature of the pref actor that will prove 
to be of some use shortly is its connection with the renormalized mean-field susceptibility, 
X- It is straightforward to verify that the second derivative of the effective potential with 
respect to the order parameter is inversely related to the susceptibility, i.e., 

n"(M) = ^ . (4.3) 

When the effective potential displays the effects of critical fiuctuations, the derivation 
above acquires modifications, specifically the inclusion of counter-terms in a full renormaliza- 
tion group treatment. Most of the steps leading to this finite size correction are as displayed 
above. The final result preserves the relationship between the prefactor and the isothermal 
susceptibility implied by the above equation. Specifically, the new prefactor has the form 



H" 1 

V 2^" V2¥^' ^^'^^ 
At the level of renormalized mean-field theory, one can write for the susceptibility 

X = e^^-")^* = e^^*/^ (4.5) 

where 7 is the critical exponent for the isothermal susceptibility (x oc |t| ''^). As noted 
previously, the quantity i* is determined by Eq. (3.3). 

There is a more general argument for the existence of the prefactor. This argument is 
based on the "infinitesimal" momentum shell version of the Renormalization Group used 
in the calculation of the partition function of the Ginzburg-Landau- Wilson model [18]. In 
this approach, modes at the surface of a shrinking Brillouin zone are integrated out. The 
contributions to the free energy have a gaussian-like form. The key contribution to the net 
free energy is 

l^log[E, + ^2] (4.6) 

k 

In expression (4.6), the self energy term contains the effects of fluctuations whose wave 
vectors exceed k in magnitude. If the uniform mode is singled out, the sum in (4.6) is carried 



out over all non-zero /c's. The one term that does not contribute to the free energy is, thus, 
equal to 



In this term the restriction on the momenta of the fluctuations contributing to Eo no longer 
differentiates the quantity from the standard self energy, E. Standard arguments suffice to 
establish the connection between the the susceptibility and the self energy. If the A; = 
term is added to the sum in Eq. (4.6), we end up with a sum over all k that constitutes the 
leading order bulk contribution to the effective potential. It is then necessary to perform a 
subtraction leading to a term in the exponent of the form 1/2 log [Eq] = Iny^xr- That the 
"subtracted" term appears in the exponent with a positive sign reflects the fact that the 
partition function is the exponential of minus the free energy. 



There are other finite size effects. For example, there is a limit on the maximum possible 
size of a block spin in the renormalization procedure. This limit appears as a modification of 
Eq. (3.3). A limit on the block spin size, and hence the quantity £*, is enforced if Eq. (3.3) 
is modified so as to ensure that the size of the block spin does not exceed the dimensions of 
the system. The new requirement on £* is 



The quantity c is a number of order unity. In principle, the results of a calculation will not 
depend on the specific value of c, as the value of £* is a detail of the renormalization group 
procedure, which has no effect on the final result [14]. In practice, the results that follow from 
the use of Eq. (5.1) depend sensitively on the choice of the parameter c. In the work reported 
here, this quantity thus acquires the status of a fitting metrical parameter connecting the 
lattice size in the Ising simulations to the length parameter L in the renormalization group 
approach. 




(4.7) 



V. FINITE SIZE EFFECTS 




(5.1) 
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VI. COMPARISON WITH DATA 



The effective potential (3.2) is now compared to the results of the Ising model simulations 
performed by Tsypin [1-3]. Four data sets from Tsypin's simulations are fit, one in the 
ordered phase {(3o — 0.2227), one in the disordered {(3d — 0.22055) and two at the critical 
temperature T — T^, {Pc — 0.22165). The data in the disordered phase are gathered in a 
Monte Carlo investigation of the order parameter distribution of a system of 58^ Ising spins. 
In the case of the ordered phase, the simulations were performed on a system of 74^ spins, 
while lattices of 16^ and 32^ were used for the simulations at criticahty. The data are fit with 
the renormalized mean-field approximation Eq. (3.2) including the prefactor. The fitting 
procedure includes the use of the restriction on i* obtained by solving the transcendental 
equation, (5.1). The overall normalization is set such that the area under our curve matches 
that of Tsypin's histograms. 

Our fitting procedure utilizes a grid-search minimization of the x^((5) merit function. 

X^(«) = ^E^^^^^^4^ (6.1) 

The parameters {aj} for the zeroth-order fit are t, u and c that appear in (3.2) and (5.1). 
Technically, this expression is the reduced since we divide by the number of degrees of 
freedom, D. D S& the number of data bins less the number of free parameters. The errors Ui 
were assumed to be of Gaussian order; i.e. cr| = y{xi\ a). This assumption is conservative 
since the is less than the expected result, which should be approximately unity. 

As an initial effort to minimize ^ blind fit is performed across all the histograms of 
a given data set (e.g., in the disordered phase the seven histograms corresponding to the 
seven different values of the external field h). See Table I. In this blind fit each bin of each 
histogram is treated with the same weight. This procedure certainly leads to an unbiased 
minimization of x^, but further fine-tuning is required as explained below. 

It is important to note the similarity between the values of u in both temperature regimes. 
In zeroth-order approximation we find ua — Ug — 0.231 for the disordered and ordered 
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phases, respectively. This is quite significant in that Landau-Ginzburg theory dictates that 
the fourth-order couphng be the same everjrwhere. Also note that the value of t that best 
fits the data from the critical simulations is indeed zero, as expected. We do observe one 
anomaly in the best fit values of the reduced temperature in the ordered and disordered 
phases, however. Namely, if we compare the ratio of the temperatures we obtain 

— = 0.873 

-to 

while from Tsypin's simulation we have 

td _ Pd — Pc _ ]^ Q5 
to Po Pc 

This is unfortunate as we expect this ratio to be the same in the simulations and our fit 
even though the temperatures themselves need not be. 

An additional comment needs to be made for the fitting of the ordered phase data. 
Notice that in the h — histogram we have removed the data lying to the left of the peak. 
Attempts to fit the data in that regime did not produce satisfactory results. We attribute 
this failure to the fact that a single renormahzed Ginzburg-Landau form (such as Eq. (3.2)) 
is inappropriate to the region of two-phase coexistence. 

One expects that the scaling form for the free energy (combined with the prefactor) 
should provide a reasonable approximation to the free energy in all temperature regimes 
using a single value of u and c. To this end, we search for a set of parameters that minimize 
the ioT all data sets. The most efi^ective strategy for accomplishing this task is to first 
determine a range of parameters u and c that fit the t — data. These bounds are then used 
for a three parameter fit in the t > and t < phases which yields final values oi t, u and 
c. It should be noted that a significant amount of "tweaking" of the parameters is necessary 
in addition to the blind minimization of the x^. These values are tabulated in II and the 
resulting probability densities are plotted against Tsypin's data in the three temperature 
regimes. The plots are displayed in Figures 1-3. It is interesting to note that the ratio of 
temperatures discussed above is now 
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1.06 



-to 

which is within one percent of Tsypin's value. 

To further evaluate our scaling form we also calculate the universal quantity [19,7] 

14 = JMXr - 3 (6.2) 

The values of r4 from our calculations compared to other sources are tabulated in III. 

As a final check on our results, we compare our best fit value of the renormalized coupling 
constant u with results obtained by other methods. A variety of techniques have been 
utilized to determine this universal quantity, including the e expansion [21], high order loop 
expansions in three dimensions [22-24] , Monte Carlo simulations and high temperature series 
[25] . Recently reported values of the renormalized fourth order coupling constant range from 
0.233 [23] to 0.236 [22]. This is to be compared with our optimum global fit m = 0.228; see 
Table II. Given that no attempt twas made to match previously determined values of the 
renormalized fourth order coupling in our fitting procedure, the quality of agreement can be 
described as, at the very least, encouraging. 
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TABLES 

TABLE I. The values of fitted parameters t, u and c obtained from the blind fits for four 

different data sets. The stated errors are 68% confidence limits; i.e. one standard deviation. 

T = Tc T = Tc T>Tc T <Tc 

L 16 32 58 74 

. n +0-0002 n +0-00006 n nnQc;«+0-0002 n nn/i i s+o ooo3 

-0.0002 U -0.00007 "-'•"-'UoOD_o.oo03 ^•^'-'^-'-"-O.OOOS 

994+0-005 917+0-005 o9'^1 +0-004 r,o-. +0.003 

U U.ZZ41:_Qoo7 ^•^-'-'-0.005 ^-^^^^-0.004 ^-^^^^-0.003 

0.437 0.458 1.1% 1.60 1?! 

0.000335 0.000237 0.0170 0.0388 

^ Due to the extremely weak dependence of the free energy on the variable c in the T phases 
there is a large range of values within the 68% confidence limit. 

TABLE II. The values of fitted parameters t, u and c after constraining to fit all temperature 
regimes simultaneously. 

T = Tc T>Tc T <Tc 

L 16, 32 ^ 58 74 

. PI +0.0001 PI nnzLns+0-0003 ^ nn'^sf^+o-ooo2 

* U_o.oo03 U.004U8_o 0006 -U.00385_o.ooo4 

n 998+0.006 99Q+0.008 pi 990+0.006 

a U.zzo_o,o23 ^-^^"-0.005 "-'■^'^°-0.003 

r n 45 n 45 +2-5 A At: +2.5 

C U.40 _o.o7 _o 45 U.40 _q^^ 

0.00112 0.0366 0.0703 
^ We fit both L = 16 and L = 32 data sets simultaneously. 
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TABLE III. The value of the universal quantity r4 calculated for c/ = 3 Ising systems at t = 



from various sources. 



source 



r4 



Independent Fit 
Independent Fit 
Constrained Fit 



16 
32 
16,32 
Simulations 



-1.4209 
-1.4074 
-1.4168 



Tsypin [3] 
Tsypin [3] 
Barber, et. al. [20] 
Barber, et. al. [20] 



16 
32 
16 
32 



-1.424(3) 
-1.410(3) 
-1.4239(6) 
-1.4095(18) 



' Here we calculate r4 using the parameters u = 0.228 and c = 0.45. 
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VII. CONCLUDING REMARKS 



We arrive at the conclusion that a low order scahng form provides an excellent fit to 
the order parameter distribution both quite near the three-dimensional Ising model critical 
point and exactly at criticality. Our approximant, while at the lowest nontrivial order in 
the interdimensional e expansion, nonetheless includes full renormalization group fiows for 
the relevant variables and thus has critical point (hypcr-)scaling built in. The scaling result 
that at criticahty H ~ |M|2<^/('^-2+'?), where 2d/{d - 2 + rj) ~ 5.7 is in close agreement with 
Tsypin's [2,1] polynomial fit. That this ought to be so follows from straightforward scahng 
arguments. 

The fitting parameters used were the nonuniversal metrical factors for the reduced tem- 
perature, t, the fourth order coupling constant u and the parameter c, which controls the 
effects of finite size on the determination of the quantity £* through the relationship Eq. 
(5.1). This is certainly fewer than allowed; the 4>^ model and the Ising model are presum- 
ably in the same universality class, but are not identical. We also took the coupling u to be 
constant. However, there is no reason to beheve that the Ising model at criticality is at its 
fixed point for coupling constants. Hence, corrections to scaling should certainly be allowed, 
which in practice would mean the renormalization group fiow for the fourth-order coupling 
u{i) could be included. We have chosen not to do so in order to minimize the number of 
free parameters. 

The prefactor used by Binder [16] and Tsypin [1,2] is discussed in Section IV within the 
setting of a renormahzation group calculation. While this prefactor has only hmited effect 
on the quality of fits in the high and low temperature phases, it plays a crucial role at the 
critical point. This issue is discussed in more depth in the Appendix, where the infiuence of 
the prefactor on the universal ratio r4 is also explored. 

The fact that a low-order-in-e expression for the renormalized free energy of the Ising 
system reproduces simulation data to such a high degree of accuracy has encouraging impli- 
cations with respect to the derivation of an equation of state applicable to an 0(1) system 
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in the immediate vicinity of a critical point. In addition, there is every reason to be hopeful 
that low order corrections to the zero loop expression utilized in this work will allow for an 
even higher precision fit to data. 
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APPENDIX A: EFFECT OF THE PREFACTOR ON r4 



The pref actor represents a finite size correction to the free energy. To see that this is so, 
it suffices to exponentiate the prefactor. The contribution of the prefactor to the free energy 
is nominally independent of the size of the system, while the leading order terms scale as 
the system's volume. Clearly, in the thermodynamic limit, the prefactor is swamped by 
those terms. Nevertheless, the prefactor cannot be ignored, even as the volume approaches 
infinity. That this is so is evident when one plots out the order parameter distribution at 
the critical point. Given Eq. (2.3), neglect of the prefactor leads to an order parameter 
distribution at the critical point {t = 0, h = 0) that has the following form 

P(M) oc e-^^'l^l^'^'™ . (Al) 

Such a distribution does not manifest one important feature of the critical point order 
parameter distribution, and that is a depression in the in the neighborhood of M = 0. The 
absence of this feature has an effect on the fit to data, and in addition, on the quantity 
defined in Eq. (6.2). As an indication of the importance of the prefactor on this universal 
ratio, we first calculate its value in the absence of this contribution. We make use of the 
general relationship 

/°° M^^e"^^'"^^""'^"^ ^ ^/[-(p+i)2d/(d-2+n) d-2 + r] r°° ^^p_^_^-)^a-2+,j)/2d-i^-y^y 
Jo 2d Jo 

^ ^-2 + ?7 ,-fa+i)2d/(d-2+r,)p ( {p+l){d-2 + r]) \ 

2d \ 2d J ' ^ ^ 

where the function F is the standard gamma function. The universal combination in Eq. (6.2) 
is, then, equal to 

r(at2M)r(^) 



Y ^ 3(d-2+rj) ^ 



3 . (A3) 



When d — 3 and r) — 0.037, one finds r4 = —0.987593. As a comparison, the unrenormahzed 

Ginzburg-Landau theory predicts r4 = —0.81156. By contrast, numerical studies indicate 
r4 ~ —1.4. [3,20] While renormalization of the theory changes r4 in the right direction, 
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the full scaling form is inadequate to the challenge of reproducing the correct value of this 
quantity. 

The defect in the above analysis lies in the fact that finite size effects have been ignored. 
These effects are crucial at the critical point, where the bulk correlation length is infinite. 
As it turns out, the most important finite size effect is the prefactor. The infiuence of 
the prefactor is dramatically highlighted when one ignores the finite size contribution to 
Eq. (5.1) for t. If the coefficent c is set equal to zero, the one finds immediately 

e^' ^ {^umY"'-'^"^' . (A4) 

This means that the order parameter distribution is of the form 

P(M) oc m(2-'')/('^-2+^) exp \-,^M^<^l^'^-^+'>^] (A5) 

The quantity k in the above equation is a combination of the size of the system, L, and the 
fourth order coupling constant, u. This quantity scales out of the result for r4 in the same 
way that the constant A divided out in the ratio in Eq. (A3). The result for r4 following 
from Eq. (A5) is then 



V 2d J \ 2d 



-3 . (A6) 



* p ^ 2-T;+3((i-2+??) j 

Inserting appropriate values for d and 77 into Eq. (A6), we obtain a r4 = —1.69147. In 
utilizing the prefactor as the sole finite size correction, we have overshot the proper value of 
r4. Note, however, that the value obtained is independent of system size. Even though the 
prefactor is a finite size correction to the free energy, in the case of r4 it is fully as imporant 
as the "leading order" contributions to the order parameter distribution. 

As for the distribution P(M), it is immediately evident that the prefactor leads to a 
depression in the neighborhood of M = 0. In fact, the order parameter distribution is forced 
to go to zero at M = by the prefactor in Eq. (A5). This is a more pronounced effect 
than is desired. The supression of the order parameter distribution is reduced when when 
finite size effects are restored to Eq. (5.1). It is the combination of these latter finite size 

20 



effects and tfie finite size correction embodied in tfie prefactor tliat yields a proper order 
parameter distribution at tfie criticai point and a vaiue of r4 tliat is in agreement witli 
previous determinations of tliis quantity. 
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FIG. 1. The order parameter distribution in the disordered phase, as obtained by Tsypin [1] 
(data points) and the results of the best constrained fit based on the effective Hamiltonian embodied 
in Eqs. (3.2) and (5.1). 
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2. The order parameter distribution in the ordered phase, as obtained by Tsypin[l] (data 



points) and the results of the best constrained fit based on the effective Hamiltonian embodied in 
Eqs. (3.2) and (5.1). 
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FIG. 3. The order parameter distribution of the system at i = 0, as obtained by Tsypin[3] and 
the results of the fit using Eqs. (3.2) and (5.1). We include both the 16^ and 32^ systems on the 
plot. 



24 



